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Abstract 

With the help of the boost operator we can model the interaction between a 
weakly interacting particle(WIMP) of dark matter(DAMA) and an atomic nuclei. 
Via this "kick" we calculate the total electronic excitation cross section of the helium 
atom. The bound spectrum of He is calculated through a diagonalization process 
with a configuration interaction (CI) wavefunction built up from Slater orbitals. All 
together 19 singly- and doubly-excited atomic sates were taken with total angular 
momenta of L=0,1 and 2. Our calculation may give a rude estimation about the 
magnitude of the total excitation cross section which could be measured in later 
scintillator experiments. The upper limit of the excitation cross section is 9.7 • 10~ 8 
barn. 
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1 Introduction 

Searching for DAMA with WIMP is an interesting question from both the- 
oretical [1] and experimental sides. A considerable experimental work are in 
progress to measure DAMA-nucleus interaction in different scintillation set- 
ups such as liquid xenon [2], Nal [3] or different anisotropic crystals [4,5]. 
More technical details about the running experiments can be found under 
[6]. Theoretical considerations state that the WIMP can scatter from a nu- 
cleus either via a scalar (spin-independent) interaction or via an axial-vector 
(spin-dependent) interaction [7]. 
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In the following we use the binary encounter approach and apply the boost 
operator, to model the "kick" between the unknown WIMP-DAMA and the 
nucleus of the helium atom and calculate the total electronic excitation cross 
section. We use a CI wavefunction built up from Slater-like orbitals to de- 
scribe the ground state and the low-lying bound spectrum of helium. Our 
CI wavefunction was successfully applied to describe different time-dependent 
problems such as heavy-ion helium collisions [8,9] or photoionization of he- 
lium in short XUV laser fields [10]. According to our knowledge, there are no 
theoretical calculations from this type forecasting measurable total excitation 
cross sections. 

Atomic units are used throughout the paper unless otherwise mentioned. 



2 Theory 



At first we have to calculate the low-lying bound spectrum of the He. We ob- 
tain the eigenfunctions and the eigenvalues by diagonalizing the time-independent 
Schrodinger equation 

H He ^ = Efr, (1) 

where Hue is the spin independent Hamiltonian of the unperturbed helium 
atom 

A ,. = M + Ei_2_i + * (2) 

2 2 ri r 2 |ri - r 2 | 
and $j is the CI wavefunction built up by a finite linear combination of sym- 
metrized products of Slater orbitals 

0(r 1 )=c(n,/ t K- 1 e-^, m (^,^), (3) 

where c(n, k) is the normalization constant. We use Slater functions with an- 
gular momentum I = 0, 1, 2 and couple them to L — 0, 1 and 2 total angular 
momentum two-electron states. In our basis we apply 9 different s orbitals, 
6 different p orbitals and 4 different d orbitals, respectively. Table I presents 
our bound He spectrum compared to other, much larger ab initio calcula- 
tions [11,12,13]. We implement the complex scaling [14] method to identify 
the double-excited states in the low-lying single continuum. It is well known 
that the lsls ground state is highly angular correlated, and further pp and dd 
terms are needed to have a accurate agreement with experimental data which 
is —2.904 a.u. [10]. We checked the role of these terms and found that the 
affect in the final total cross sections is negligible. 

We may approximate the interaction between the unknown DAMA particle 
and the nucleus of the helium with the boost operator. If we suddenly "kick" 
the He nucleus with a k boost in the direction of r that is equivalent to a 
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collective — k/2 "kick" of the two atomic electrons according to the center-of- 
mass. For a better understanding the geometry of the interaction is presented 
in Figure 1. 

The total excitation amplitude can be calculated in the following way: 

a exc = ^($(r 1 ,r 2 ) / |e- irik / 2 - ir2k / 2 |$(r 1 ,r 2 ) lsls ), (4) 
/ 

where lsls is the ground state of He and the summation / runs over the 
singly- and doubly-excited final states. For elastic collision only the ground 
state to ground state transition is considered. The energy of the unknown 
DAMA is E = k 2 /2. The DAMA-electron interaction is evaluated thought the 
transition matrix elements of the boost operator between two Slater orbitals 
(0i(r|e _o " k//2 |0 2 (r)). To separate the radial and the angular part of the matrix 
element we expand the plane wave through spherical Bessel functions in the 
well-known way [15]: 

oo +1 

e irk = e k - rcos ^ =AttJ2 E i l 3i{kr)Y* m {e^^)YiA0r^r)- (5) 

1=0 m=-l 

After some algebra the angular part of the matrix element gives us the Clebsch- 
Gordan coefficient 



jY* imi (9 : <p)Y hm (9,<p)Y l2:m2 (9,<p)dn 



x 



(2/ 2 + 1)(2/ + 1) 
\ 2ir(2h + 1) 

(Z 2 , 1, h\m 2 , m, mi) ■ (Z 2 , /, /i|0, 0, 0) 



(6) 



According to the definition of the spherical Bessel functions [17] a ji{kr) = 
\f^Ji+i/2{kr) the radial part of the matrix element has the analytic solution 
of [16]: 



oo 

J J v (kr)e- ar r' t - 1 dr = 



x 



2^1 



r(v + i)^/(k 2 + a 2 y+» 

u + fil— fi + u _ k 2 



k 2 + a 2 



(7) 



where T is the gamma function and is the hypergeometric function with the 
following real arguments fi — 1 = n 1 +n 2 — 1/2, v = 1 + 1/2 and a = Ki+k 2 , We 
tried to simplify the final formula but unfortunately, we could not succeed, 
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and have to calculate the hypergeometric function numerically with a well- 
behaving complex contour integral [17]b. It is worth to mention that with 
additional constraints among the parameters [a, k, v, /i] this radial integral 
can be simplified, but not in our general case. 

The total excitation cross section can be evaluated with the following formula 

(Texc = r 2 a lxP exc , (8) 

where r a = 1.76 • 10~ 15 m is the radius of the He nucleus and P exc = \a exc \ 2 is 
the total excitation probability. For elastic collision only the ground state to 
ground state transition is considered. 



3 Results and discussion 

Figure 2 presents our elastic and excitation total cross sections in function of 
the impulse of the DAMA. The cross sections are given in barns and the wave 
number of the unknown particle is given in atomic units. If the velocity of the 
DAMA is known then the mass can be calculated from m = k/v. 

The maximum of the excitation cross section is 9.7 x 1CT 8 barn at k=3 a.u. 
DAMA impulse. At low wave numbers (k) the boost operator can be well 
approximated with its Taylor series which is similar to the dipole interaction, 
and used in photoionization calculations. At low k values the elastic cross 
section is many magnitude higher than the excitation one, which meets our 
physical intuition. At larger wave numbers, however, the dipole approximation 
breaks down and the general matrix element have to be calculated where the 
non-dipole contributions play a significant role. Above k = 30 a.u. the elastic 
and excitation cross sections run together, but the excitation cross sections 
are a factor of 2-4 higher than the elastic ones. At wave numbers larger than 
10, due to the quick oscillations of the boost operator, the cross sections have 
a strong decay which can be excellently fitted with the following power law: 

cw = 2.3566 x 10 5 • AT 13 ' 782 , (9) 

where the standard error of the exponent is 0.074 and the standard error of 
the scaling constant is 0.48, respectively. 

At k > 10000 wave numbers the cross sections stop decaying and show spurious 
oscillations, which are numerical art-effects due to the limited accuracy of the 
calculations. These cross sections are not presented in Fig. 2. 

We can not enhance the total angular momenta of the two electron wavefunc- 
tion in our calculation, and the number of the available bound states are also 
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limited. The role of the highly-excited Rydberg states are out of our scope 
too. In this sense we can not rigorously prove the convergence of our calcula- 
tion, but our experience shows that for excitation the significant contributions 
always came from the lowest excited states. We interpret our results as a rude 
approximation for the dark matter He interaction which may stimulate further 
investigations. 



4 Summary and Outlook 



With the help of the boost operator we gave a "simple-man's model" for the 
DAMA-helium nucleus interaction and calculated total electronic excitation 
and elastic collision cross sections which can be measured in future scintillator 
experiments. Our calculation could be generalized for atoms with electrons 
more than two, (even for Xe) if the wavefunction of the ground state and a 
significant large number of excited states are present with sufficient accuracy. 
We think that this problem could be solved with the General Relativistic 
Atomic Structure(GRASP) code [18], which is out of our capability. The aim 
of this paper was twofold. First, we presented our model for the DAMA- 
He interaction calculating cross sections. Secondly, we advised our model to 
many-electron-atom theorists to calculate DAMA-Xe interaction. 
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Fig. 1. The geometry of the DAMA-He interaction. 



Table 1 

The energy levels of bound, singly- and doubly-excited states used in our calcula- 
tions, compared with a) basis set calculations from [11] b) CI calculation results 
[12] and c) complex-coordinate rotation calculations from [13]. 



states 


our 


other 


states 
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other 


states 


our 
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results 


theory 




results 


theory 




results 


theory 


lsls 


-2.8821 


-2.9037" 


ls2p 


-2.1233 


-2.1238 6 


ls3d 


-2.0556 


-2.0556 6 


ls2s 


-2.1441 


-2.1460 a 


ls3p 


-2.0550 


-2.0551 6 


ls4d 


-2.0312 


-2.0313 6 


ls3s 


-2.0607 


-2.0612" 


ls4p 


-2.0259 


-2.0310 6 


2s3d 


-0.5597 


-0.5692 c 


ls4s 


-2.0333 


-2.0335 a 


2s2p 


-0.6572 


-0.6931 c 


2s4d 


-0.5305 


-0.5564 c 


2s2s 


-0.7297 


-0.7779" 


2s3p 


-0.5821 
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2s3s 


-0.5711 
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-0.5401 
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-0.5449" 


2s5p 


-0.5225 


-0.5470 c 








2s5s 


-0.5133 


-0.5267" 


3s3p 


-0.2998 


-0.3356 c 
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Impuls of the WIMP DAMA in a.u. 

Fig. 2. Elastic (dashed line) and excitation (solid line) total cross sections for 
DAMA-He collision. 
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